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Abstract 

The  ability  to  accurately  predict  the  drag  forces  on  a  ship  before  it  is  built  would  lead  to  more 
efficient  designs.  To  do  this,  the  effect  of  surface  roughness  on  frictional  drag  must  be  well 
understood.  The  goal  of  this  project  is  to  identify  the  appropriate  roughness  scaling  parameters 
for  simple  three-dimensional  roughness  with  similar  length  scales  to  those  found  on  ship  hulls. 
In  order  to  study  and  ultimately  predict  the  effects  of  surface  roughness  on  fluid  flow  and  drag, 
flat  plates  with  smooth  and  rough  surface  conditions  were  tested.  Mesh  and  sandpaper  were 
chosen  as  the  rough  surfaces  because  they  are  three  dimensional,  which  is  characteristic  of 
naturally  occurring  surfaces.  The  first  phase  of  testing  involved  towing  the  plates  in  a  tow  tank 
to  determine  the  overall  frictional  drag.  These  tests  were  done  in  the  115  m  long  tow  tank 
located  in  the  USNA  Hydromechanics  Laboratory.  Detailed  velocity  measurements  were  also 
obtained  with  the  plates  in  a  re-circulating  water  channel  located  in  the  Hydromechanics 
Laboratory  to  determine  the  effect  of  the  roughness  on  the  turbulence  near  the  surface.  The 
velocity  measurements  were  obtained  with  a  laser  Doppler  velocimeter  (LDV).  The  drag  results 
from  both  sets  of  tests  showed  excellent  agreement.  It  was  also  observed  that  beyond  a  few 
roughness  heights  from  the  wall,  the  normalized  turbulence  was  independent  of  the  roughness. 
Proper  scaling  for  the  sandpaper  was  found  to  be  a  function  of  the  roughness  height,  while  the 
mesh  surfaces  were  discovered  to  be  a  function  of  both  roughness  height  and  wire  spacing.  The 
results  from  this  project  will  aid  in  the  development  of  a  general  model  of  overall  frictional  drag 
from  physical  measures  of  the  surface  alone. 
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1.  Nomenclature 
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2.  Introduction 

Being  able  to  aeeurately  prediet  the  foree  of  drag  on  a  ship  or  an  airplane  before  it  is  built 
would  greatly  aid  in  the  design.  The  United  States  Navy  would  benefit  greatly  from  this  ability. 
Predieting  the  drag  forees  on  a  ship  before  it  is  built  would  allow  for  them  to  be  designed  to  sail 
faster,  longer,  and  more  effieiently.  Furthermore,  understanding  and  predieting  the  effeets  of 
drag  would  help  in  the  development  and  performanee  of  turbines,  eompressors  and  other  bladed 
turbomaehines  whose  eharaeteristies  are  ehanged  and  affeeted  by  surfaee  roughness(Aeharya, 
Bornstein,  and  Eseudier  1986). 

When  designing  ships,  one  method  of  predieting  drag  involves  model  testing.  The  laws 
of  similitude  allow  for  the  behavior  of  full  seale  prototypes  to  be  predieted  from  models.  For 
model  ships,  it  is  important  that  the  model  and  prototype  have  the  same  shape  and  differ  only  in 
size.  One  problem  is  that  the  size  of  the  surfaee  roughness  is  very  hard  to  reduee  proportionally 
to  the  redueed  size  of  the  ship  (Franzini  1997).  Henee,  it  is  impossible  to  prediet  how  mueh  of 
the  drag  is  aetually  eaused  by  the  surfaee  roughness.  In  order  to  eompensate  for  this,  engineers 
add  a  eorrelation  allowanee,  Ca,  to  the  total  resistanee  eoeffieient,  Cr,  to  get  a  drag  eoeffieient 
whieh  aeeommodates  surfaee  roughness  (Fewis  1988).  The  allowanee  eoeffieient  is  eompletely 
empirieal,  representing  an  average  roughness  eondition.  Further  work  was  done  speeifieally  in 
hull  roughness  effeets  whieh  related  a  roughness  height  to  Ca  (Bowden  &  Davison  1974).  The 
problem  with  this  method  is  that  the  effeetive  roughness  height  still  needs  to  be  found 
experimentally.  An  improved  resistanee  eoeffieient  would  relate  surfaee  frietional  drag  to  a 
speeifie,  measured  surfaee  roughness.  This  knowledge  would  greatly  aid  ship  designers. 

The  purpose  of  this  researeh  projeet  was  to  examine  the  effeets  that  two  different  surfaee 
roughness  types  have  on  drag.  From  the  results,  a  model  whieh  eorrelates  these  partieular 
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surface  roughnesses  to  their  assoeiated  drag  ean  be  created.  The  model  results  eould  then  be 
scaled  up  to  ship  seale  to  see  how  the  surface  roughness  inereases  the  overall  drag  and  required 
propulsion  power  of  a  ship. 
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3.  Background 

Drag  is  defined  as  the  net  foree  in  the  direetion  of  the  flow  opposing  the  motion  of  the 
body  through  a  fluid  (Alexandrou  2001).  For  a  ship,  drag  is  eomprised  of  two  main  eomponents, 
wave-making  drag  and  viseous  drag,  whieh  ean  be  added  together  to  get  the  total  drag  foree. 
When  ships  move  through  water,  waves  are  ereated.  These  waves  are  ereated  by  the  varying 
pressure  that  the  ship  is  eausing  under  the  water  and  the  eonstant  pressure  at  the  surfaee  (Tupper 
1996).  The  energy  needed  to  make  these  waves  is  provided  by  the  ship.  Viseous  drag  is 
eomprised  of  both  pressure  drag  and  frietion  drag.  Pressure  drag  is  the  drag  produeed  by  normal 
stresses  on  the  surfaee  of  the  body.  A  plate  perpendieular  to  the  flow  has  nearly  all  pressure  drag 
whieh  is  eaused  by  the  differenee  in  pressure  between  the  flow  on  the  upstream  and  downstream 
side  of  the  plate  (Munson,  Young,  and  Okiishi  2002).  For  a  plate  parallel  to  the  flow,  the 
pressure  drag  is  minimal,  and  the  drag  is  mainly  due  to  frietion.  Frietion  drag  is  ereated  by  shear 
stresses,  whieh  are  eaused  by  viseous  and  turbulent  effeets.  With  turbulent  flows,  whieh  inelude 
most  flows  of  praetieal  interest,  surfaee  roughness  ean  signifieantly  inerease  the  frietional  drag 
foree  (Munson,  Young,  and  Okiishi  2002). 

Numerous  researeh  studies  (Sehoenherr  1932,  Granville  1987,  and  Grigson  1987,  1992) 
have  been  eondueted  regarding  how  surfaee  roughness  affeets  a  ship’s  drag.  These  studies  foeus 
on  ways  to  seale  tow  test  data  of  models  eovered  with  a  surfaee  roughness  up  to  that  of  full  seale 
ships  in  order  to  determine  the  aetual  drag  assoeiated  with  surfaee  roughness.  Beeause  model 
tests  alone  will  not  give  the  full  seale  viseous  drag  on  a  ship,  these  studies  have  foeused  on 
measuring  the  assoeiated  roughness  funetion  whieh  ean  then  be  used  to  get  the  full  seale  viseous 
drag  on  a  ship  (Grigson  1987).  Other  methods  for  predieting  the  behavior  of  Aids  on  surfaces 
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involve  using  Computational  Fluid  Dynamics  (CFD)  to  perform  Direct  Numerical  Simulation 
(DNS).  The  problem  with  DNS  in  the  present  case  is  that  computing  turbulent  flows  requires  a 
very  large  computational  grid  at  high  Reynolds  numbers.  There  are  two  main  length  scales  that 
DNS  have  to  account  for;  the  viscous  length  scale  which  is  of  the  order  of  microns,  and  the 
boundary  layer  length  scale  which  is  of  the  order  of  centimeters.  The  DNS  must  capture  the  flow 
physics  at  the  viscous  length  scale  throughout  the  entire  boundary  layer.  At  present,  computers 
are  not  able  to  perform  the  large  number  computations  in  a  reasonable  time.  Some  CFD  codes 
using  turbulence  models  using  a  simplified  computational  grid  including  numerous  assumptions 
have  been  developed  which  includes  surface  roughness  effects.  Furthermore,  results  from 
experiments  are  still  needed  as  inputs  to  the  models. 

As  fluid  flows  over  a  surface,  a  thin  layer  develops  right  above  the  surface  due  to 
viscosity.  In  this  thin  layer,  called  the  boundary  layer,  the  velocity  of  the  flow  varies  from  zero 
at  the  wall  to  the  flow  speed  at  the  edge  of  the  layer.  An  example  of  a  turbulent  boundary  layer 
mean  velocity  profile  can  be  seen  in  Figure  1 .  Many  things  affect  this  thin  boundary  layer, 
including  the  type  of  flow  (laminar  or  turbulent),  the  speed  of  the  flow,  and  fluid  properties. 

With  turbulent  flows,  surface  roughness  affects  the  boundary  layer  near  the  wall  by  creating 
higher  wall  shear  stress.  This  in  turn  creates  more  frictional  drag. 
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u/u, 
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Figure  1-Turbulent  boundary  layer  velocity  profile:  (a)  boundary  layer  development; 

(b)  non-dimensional  profile 


A  turbulent  boundary  layer  ean  be  modeled  as  a  composite  of  various  layers,  as  shown  in 
Figure  2.  The  first  layer  adjacent  to  the  wall,  which  is  known  as  the  viscous  sublayer,  has  shear 
stress  caused  almost  entirely  by  viscosity.  Viscosity  is  a  fluid’s  internal  resistance  to  shear 
deformation  (Alexandrou  2001).  This  sublayer  is  extremely  thin,  generally  extending  only  a  few 
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hundredths  of  a  millimeter  out  from  the  wall  (Franzini  1997).  The  outer  layer  is  the  region  in 
whieh  viseosity  does  not  have  a  direet  effect  on  shear  stress.  Instead,  the  shear  stress  in  this 
outer  layer  region  is  caused  by  turbulent  motions.  Between  the  sublayer  and  outer  layer,  the  log- 
law  region  is  the  location  where  the  total  shear  stress  is  nominally  constant  and  maximum. 
Turbulent  stresses  also  dominate  in  this  region.  The  log  region  is  characterized  by  the  linear 


velocity  profile  when  plotted  in  log-normal  coordinates  (Clauser  1954). 


Figure  2-The  mean  velocity  profile  for  a  smooth-wall  turbulent  boundary  layer 


Typically,  a  boundary  layer  over  a  rough  surface  is  modeled  by  the  same  methods  as  a 


smooth  wall  boundary  layer.  This  model  uses  a  standard  two-layer  approach  in  which  only  the 


inner  layer  is  directly  affected  by  the  roughness,  while  the  outer  layer  is  affected  only  indirectly 


through  the  increase  in  the  wall  shear  stress.  The  effects  of  surface  roughness  close  to  the  wall 
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are  not  well  understood  (Aeharya,  Bornstein,  and  Eseudier  1986).  Further  researeh  is  needed  in 
this  area. 

Surfaee  roughness  is  thought  to  affeet  the  mean  veloeity  profile  by  shifting  it  downward 
when  plotted  in  log-normal  eoordinates.  The  drag  assoeiated  with  the  roughness  deereases  the 
momentum  of  the  flow.  This  loss  of  momentum  ean  be  seen  in  the  Figure  3  where  Alf,  the 
roughness  funetion,  represents  the  effeet  of  roughness  on  the  mean  veloeity  profile.  The  larger 
the  roughness  funetion  is,  the  more  effeet  the  roughness  has  on  the  mean  profile  and  henee  the 
wall  shear  stress  (Hama  1954).  This  inerease  in  wall  shear  stress  will  in  turn  ereate  more  drag. 
Another  important  feature  to  note  is  how  the  rough  surfaee  profile  maintains  the  same  shape  as 
the  smooth  wall  profile.  If  the  surfaee  roughness  were  to  affeet  the  outer  layer,  the  mean 
veloeity  profile  for  the  rough  surfaee  in  the  outer  layer  would  not  have  the  same  shape  as  the 
smooth  wall.  Sinee  the  eurrent  two  layer  model  of  the  traditional  boundary  layer  assumes  that 
the  outer  layer  is  not  affeeted  by  surfaee  roughness,  modifieations  of  the  eurrent  model  would 


need  to  be  made. 
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Figure  3-Turbulent  boundary  layer  profiles:  smooth  and  rough  wall 


The  effects  of  surface  roughness  on  mean  velocity  profiles  have  been  the  topic  of  many 
research  studies.  This  can  be  seen  in  the  works  of  Clauser  (1954),  Hama  (1954),  Granville 
(1987),  and  Schultz  and  Myers  (2003).  Much  of  this  work  was  done  in  order  to  determine  the 
roughness  function  and  to  show  that  increasing  the  surface  roughness  merely  increases  the 
downward  shift  of  the  mean  velocity  profde.  For  these  studies,  the  roughness  function,  Af/^, 
was  obtained  by  measuring  the  mean  velocity  profile  directly  or  by  using  a  tow  tank  and  backing 
out  the  roughness  function  via  boundary  layer  similarity  models. 

Previous  research  has  shown  that  at  sufficiently  high  Reynolds  number,  surface 
roughness  does  not  affect  the  turbulent  flow  outside  a  sublayer  that  extends  a  few  roughness 
heights  above  the  surface  (Rotta  1962,  Townsend  1976).  This  implies  that  very  little 
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communication  occurs  between  the  near-wall  region  and  the  rest  of  the  flow.  Reeent  papers, 
however,  elaim  that  the  effeets  of  surfaee  roughness  are  not  eonfined  to  the  wall  region  (Antonia 
and  Krogstad  2001),  but  ean  extend  beyond  the  log  law  region.  Additionally,  they  elaim  that 
surfaee  roughness  affeets  the  turbulent  stresses  that  eontribute  to  surfaee  drag.  However,  there  is 
not  a  general  eonsensus  regarding  the  extent  and  influenee  of  surfaee  roughness  on  turbulenee 
among  researehers  in  the  field.  (Krogstad  and  Antonia  1999). 

Beeause  the  traditional  two  layer  model  of  the  turbulent  boundary  layer  has  reeently  been 
ealled  into  question,  new  researeh  must  be  done  in  order  to  explain  the  confiieting  findings 
whieh  have  been  seen  by  some  engineers.  Some  of  the  reeent  findings  may  have  been  a  result  of 
using  different  surfaee  roughness  geometries  whieh  affeeted  the  boundary  layer  in  a  unique  way 
that  is  not  oharaeteristie  of  other  types  of  surfaee  roughness.  Furthermore,  there  must  be  some 
roughness  height  at  whieh  eurrent  boundary  layer  models  fail  beeause  the  surfaee  roughness 
beeomes  too  large  for  its  effeets  on  the  flow  strueture  to  be  eonfined  to  the  inner  layer.  This 
information  would  be  very  helpful  in  predietions  of  drag  for  ships. 

The  goal  of  this  projeet  was  to  experimentally  investigate  these  questions.  The  drag  on 
three  types  of  surfaees  was  obtained.  A  smooth  plate  was  used  as  a  eontrol  surfaee.  Mesh 
surfaees  along  with  sandpaper  surfaees,  three  of  eaeh,  were  tested  as  the  rough  surfaees.  Mesh 
and  sandpaper  were  ehosen  as  the  rough  surfaees  beeause  they  are  three  dimensional,  whieh  is 
eharaeteristie  of  naturally  oeeurring  surfaees,  they  are  easily  repeatable,  and  there  is  some  data 
already  available  in  literature  for  these  surfaees.  The  height  of  the  roughness  elements  for  the 
mesh  and  sandpaper  were  varied  to  determine  if  typieal  boundary  layer  models  are  valid  for  large 
roughness  heights.  The  first  phase  of  testing  involved  towing  the  plates  in  a  tow  tank  to 
determine  the  frietional  drag  on  the  plate.  These  tests  were  done  in  the  115  m  long  tow  tank 
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located  in  the  USNA  Hydrodynamies  Laboratory.  Next,  detailed  veloeity  measurements  were 
obtained  with  the  plates  in  a  re-eireulating  water  ehannel  loeated  in  the  hydrodynamies  lab  to 
determine  the  effeet  of  the  roughness  on  the  turbulenee  near  the  surfaee,  by  measuring  mean  and 
turbulent  veloeities  in  the  boundary  layer.  These  detailed  veloeity  measurements  were  obtained 
with  a  laser  Doppler  veloeimeter  (LDV). 
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4.  Experimental  Procedures-Tow  Tank  Tests 

All  drag  measurements  made  on  the  plates  were  done  in  the  115  meter  tow  tank  loeated 

in  the  Hydromechanics  Laboratory  located  at  the  United  States  Naval  Academy.  The  width  and 
depth  of  the  tank  are  7.9  and  4.9  meters,  respectively.  The  towing  carriage,  which  is  powered  by 
two  300  kW  motors,  can  vary  its  velocity  between  0  and  7.62  m/sec.  With  this  research  project, 
the  speed  range  was  1.8  and  3.7  m/s.  The  plates  were  1.52  meters  long  and  0.56  meters  wide, 
with  a  thickness  of  3.8  mm.  They  were  covered  on  both  sides  with  the  given  surfaces  which 
were  held  on  using  epoxy.  The  plates  were  then  secured  to  the  towing  carriage  using  the 
apparatus  shown  in  Figure  4. 


TOWING  CARRIAGE  DIRECTION 
CARRIAGE  MOUNT 


Figure  4-Schematic  of  the  flat  plate  drag  apparatus  for  towing  tank  tests 
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A  stiffener  was  put  on  the  top  of  the  plates  so  that  the  plate  would  not  bend  when  being 
towed  through  the  water.  Two  foree  gages  were  then  attaehed  to  the  plate.  The  first  gage 
measured  the  foree  parallel  to  the  plate  whieh  was  the  aetual  viseous  drag  foree.  The  seeond 
gage  measured  the  forees  aeting  perpendieular  to  the  plate  whieh  were  the  side  forees.  These 
forees  are  not  important  to  the  drag  ealeulations;  however,  attempts  were  made  to  minimize  the 
side  forees  to  ensure  the  plate  was  parallel  to  the  flow.  This  was  done  by  towing  a  plate  several 
times  while  monitoring  the  side  forees,  and  then  making  adjustments  on  the  alignment  of  the 
plate. 

The  foree  gages  start  off  as  a  solid  bloek  of  magnetie  stainless  steel.  They  are  then 
hollowed  out  in  sueh  a  way  so  that  they  only  defleet  in  one  direetion  when  a  foree  is  applied.  A 
displaeement  transdueer  is  then  used  to  measure  the  distanee  the  gage  defleets.  The  transdueer 
has  a  fixed  eoil  and  a  ferrous  eore  that  is  free  to  move  in  the  direetion  of  the  displaeement.  A 
signal  eonditioning  unit,  (SCU),  sends  a  high  frequeney  AC  eurrent  to  the  eoil.  As  the  gage  is 
displaeed  by  a  foree,  the  ferrous  eore  moves  inside  the  eoil  and  ehanges  the  AC  eurrent.  This 
signal  is  then  sent  baek  to  the  SCU  where  it  is  eonverted  to  a  DC  voltage.  This  voltage  is  now 
proportional  to  the  displacement  of  the  gage. 

Before  setting  the  gages  on  the  plate,  they  have  to  be  calibrated  to  determine  what  DC 
voltage  is  output  when  known  forces  deflect  the  gages.  This  is  done  by  hanging  weights  off  the 
gages  and  recording  the  voltage  of  the  SCU.  Using  this  information,  a  linear  calibration  plot  can 
be  made  that  relates  force  and  voltage.  The  slope  of  this  plot  is  used  to  convert  the  force  gage 
output  voltages  into  force  measurements.  An  example  of  the  calibration  plot  can  be  seen  in 
Figure  5. 
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Figure  5-Force  gage  calibration  curve 

The  analog  signal  from  the  SCU  is  then  sent  to  an  analog  to  digital,  A  to  D,  eonverter.  After 
going  through  the  A  to  D  eonverter,  the  signal  is  then  sent  to  a  eomputer  for  data  aequisition. 

The  voltage  from  the  SCU  is  multiplied  by  the  ealibration  eonstant  gotten  from  the  ealibration 
plot  to  determine  the  measured  drag  foree.  A  bloek  diagram  of  the  different  parts  of  this  proeess 
is  shown  in  Figure  6. 


Figure  6-Block  diagram  for  force  measurements 
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5.  Experimental  Analysis  of  the  Tow  Tank  Results 

One  smooth  plate  and  six  plates  with  different  roughness  (three  sandpaper  and  three  mesh) 

were  tested.  The  drag  foree  for  eaeh  surfaee  roughness  was  measured.  Initially  the  wave  making 
forees  assoeiated  with  the  towed  plates  were  measured.  This  was  done  by  towing  a  smooth  plate 
with  2.54  em  of  it  submerged  in  the  water.  The  foree  measured  was  then  assumed  to  be  eaused 
primarily  by  wave-making  forees  although,  the  frietional  forees  were  aeeounted  for  in  the  final 
equation  by  subtraeting  this  wetted  surfaee  area  from  the  total  wetted  surfaee  area.  The  plates 
were  then  submerged  56  em  below  the  water  surfaee.  In  previous  researeh  (Sehultz,  private 
eommunieation),  wave  profdes  were  observed,  photographed,  and  analyzed  for  plates  that  were 
fully  submerged  and  for  plates  that  were  partially  submerged  2.54  em.  It  was  eoneluded  from 
these  tests  that  for  a  flat  plate,  waves  were  eaused  by  surfaee  piereing  rather  than  by  how  mueh 
of  the  plate  was  below  the  water  line.  When  towed  at  seven  different  speeds,  foree 
measurements,  whieh  were  ealled  the  raw  forees,  were  reeorded  by  the  foree  gauges  on  the  plate. 
The  previously  measured  wave-making  forees,  ealled  tare  forees,  Ftare,  were  then  subtraeted  from 
the  measured  raw  forees,  Fraw  The  result  is  the  frietional  drag  foree,  Fd,  on  the  plate,  as 
demonstrated  in  equation  1 . 


Fd  =  Fraw  -Ft 


tare 


(1) 


By  knowing  the  veloeity  that  the  plates  were  towed  at  and  the  drag  forees  that  were  measured, 
the  eoeffieient  of  frietional  drag,  Cp,  can  be  ealeulated.  This  is  done  by  dividing  the  drag  foree 


1  , 

by  the  dynamie  pressure,  ( —  pU  )  where  p  is  the  density  of  the  fluid  and  U  is  the  towing 


veloeity,  and  the  wetted  surfaee  area  of  the  plate,  S.  This  ean  be  seen  in  equation  2. 
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Cf^ 


Fq 

-pU^S 

2 


(2) 


For  a  smooth  wall,  Cf  is  a  function  of  the  Reynolds  number,  Rei.  The  Reynolds  number  is  a 
non-dimensional  number  which  represents  the  ratio  of  inertial  forees  to  viseous  forees,  and  as  the 
Reynolds  number  increases  beyond  the  critical  value  for  transition,  turbulence  increases 
(Franzini  1997).  The  Reynolds  number  is  expressed  in  equation  3 


where  L  is  the  length  of  the  plate  and  p  is  the  dynamie  viseosity  of  the  fluid. 


For  the  smooth  wall,  the  eoefficient  of  frietional  drag  decreases  with  inereasing  Reynolds 
number.  The  reason  for  this  is  that  both  viscous  and  inertia  forces  are  still  present,  but  viscous 
forces  are  deereasing.  For  a  fully  rough  wall,  Cf  is  independent  of  the  Reynolds  number  at  high 
Reynolds  numbers  because  turbulent  stresses  dominate  the  viscous  stresses,  as  illustrated  in 
Figure  7. 

When  Cfis  known,  Granville’s  Similarity  Theory  can  be  used  to  calculate  the  roughness 
funetion,  Af/^  (Granville  1987)  .  The  roughness  function  is  a  crucial  value  beeause  it  allows  for 
drag  to  be  scaled  up  to  larger  objects.  Granville’s  theory  states  that,  for  the  same  value  ofReiCF, 


2  2 

the  difference  between  - of  the  smooth  wall  and  - of  the  rough  wall  approximately 


equals  the  roughness  function,  Af/^  when  plotted  against  log(7?eiCf)  .  Similarity  theory  states 
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that  this  calculation  yields  the  same  result  as  measuring  the  downward  shift  in  the  mean  veloeity 
profile  at  the  trailing  edge  of  the  towed  plate.  An  example  of  this  is  shown  in  Figure  8. 


Figure  7-Coefficient  of  frictional  drag:  smooth  and  rough  walls 


Figure  8-Relationship  between  Alt  and  coefficient  of  frictional  drag 
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The  Karman  Schoenherr  line  whieh  agrees  closely  (within  1%)  with  the  present  smooth  wall 


results  is  used  to  get  —  of  the  smooth  wall  (Schoenherr  1932.  The  Karman  Schoenherr  line 


is  given  in  equation  4. 


0.242 


=  log(Rei  Cf) 


(4) 


Once  —  of  the  smooth  wall  is  known,  the  roughness  function  for  the  towed  plate  data  can  be 
A  Cc- 


found.  This  is  shown  in  equation  5. 


AU^  = 


JS 


Jr 


-19.7 


C. 


J  s 


C. 


Jr 


- AU^ 

K 


c. 


(5) 


Jr 


Equation  5  shows  Granville’s  similarity  law  where  k  is  the  von  Karman  constant,  which  is  equal 
to  0.41,  and  Alf"  is  the  roughness  function  slope.  The  roughness  function  slope  is  the  slope  of 
Alf"  as  a  function  of  ln(A:^).  The  roughness  Reynolds  number,  ,  is  the  ratio  of  the  roughness 
length  scale,  k,  to  the  viscous  length  scale,  v/Ur.  The  roughness  function  slope  is  also  needed  in 
order  to  obtain  k^.  This  can  be  seen  in  equation  6  which  is  also  a  part  of  Granville’s  similarity 


law. 
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D  r  V 


\Lj\ 


KlC, 


1-- 

K 


c. 


Jr 


+  — 

K 


- MJ^ 

Ik 


C. 


yv 


Jr 


(6) 


An  iterative  solution  is  required  to  solve  for  equations  (5)  and  (6)  since  A 17^  is  unknown  in  both 
cases.  To  begin  with,  Af/^  is  taken  to  equal  zero  and  the  equations  are  solved  and  Af/^  is  plotted 
versus  ln(A:^).  The  slope  of  this  line  is  found,  which  can  then  be  substituted  back  in  for  Af/^  . 
This  process  is  repeated  until  the  slope  of  the  line  and  Af/^  used  in  the  equations  converge. 

Since  it  is  impossible  to  determine  the  effect  that  surface  roughness  has  on  a  ship  through 
model  tow  tests,  Af/^  is  found  by  towing  flat  plates  which  can  then  be  scaled  up  to  determine 
what  Cf  is  for  full  scale  ships.  This  is  done  by  using  Granville’s  similarity  law  for  scaling  up  Cp 
(Schultz  1998).  Initially,  the  smooth  coefficient  of  friction,  Cpsmooth,  is  determined  for  a  ship 
using  equation  (4).  Crsmooth  is  then  plotted  against  log  (Rei).  A  new  curve  representing  Cpmugh 
is  created  by  shifting  the  Crsmooth  curve  by  a  distance  of  Af/^ /(ln(10)/K:)  in  the  positive  logiRei) 
direction  (Schultz  1998).  This  can  be  seen  in  Figure  9.  The  location  of  the  experiment  value  of 
Cfrough  from  the  tow  tank  tests  on  the  displaced  curve  is  also  shown  on  Figure  9.  As  expected, 
this  value  falls  on  the  Cpmugh  curve. 
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Figure  9-CFsmooth  and  experimental  Cp 

The  ratio  of  the  plots  length  to  the  viscous  length  scale  is  called  L^.  A  line  of  constant 
passing  through  the  lab  scale  data  is  created  using  equation  7. 


r 


ic 


This  can  be  seen  in  Figure  10. 
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Figure  10-Line  of  constant 

This  line  of  constant  is  then  shifted  in  the  positive  log  {Rei)  direction  by  a  distance  of 
\og{Lship  scale! Liab  scale)  (Schultz  1998).  Thc  intcrscction  of  the  shifted  line  and  the  Cpmugh  line 
gives  the  scaled  up  result  of  Cprough  for  a  full  scale  ship.  This  can  be  seen  in  Figure  11. 
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Figure  11-Schematic  of  Granville  scale-up  procedure 


Furthermore,  the  Reynolds  number  where  this  interseetion  oeeurs  ean  be  reeorded  from 
Figure  1 1 .  Using  the  Reynolds  number,  the  speed  at  whieh  the  full  scale  ship  would  be  cruising 
at  can  be  calculated.  This  is  done  using  equation  8 

Re,  V 


U  =  - 


L 


(8) 


where  U  is  the  ship’s  speed  and  L  is  the  length  of  the  ship  which  was  taken  to  be  150  meters. 


This  length  is  the  same  length  as  a  DDG-5 1  destroyer,  so  this  would  be  the  speed  that  this 
destroyer  was  operating  at.  Next  the  Cpsmooth  is  determined  for  a  full  scale  ship.  This  is  done  by 
using  equation  5  and  plugging  the  new  Rei  values  that  were  found  for  a  full  scale  ship.  With  this 


known,  the  percent  increase  in  the  coefficient  of  frictional  drag  due  to  roughness  can  be  found 
using  equation  9. 
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C 


Frough 


C 


Fsmooth 


c 

^  Fsmooth 


=  %  increase  in  frictional  drag  due  to  roughness 


(9) 


Ships  spend  most  of  their  time  operating  at  or  near  their  cruising  speed.  It  is  important  to 
note  how  roughness  affects  the  drag  at  this  speed.  For  a  DDG-51  class  destroyer,  this  speed  is 
around  15  knots.  With  this  information,  it  is  then  possible  to  find  how  the  percent  increase  in  the 
coefficient  of  frictional  drag  due  to  roughness  affects  the  total  coefficient  of  drag.  Using  data 
already  taken  from  model  tests  of  a  DDG-51,  the  coefficient  of  drag  due  to  friction  is  found  at 
certain  speeds  as  well  as  the  total  drag  coefficient.  Using  equation  10 

^  =%  of  drag  due  to  friction  (10) 


where  Cf  is  the  coefficient  of  drag  due  to  friction  and  Ct  is  the  coefficient  of  the  total  drag,  the 
percent  of  drag  due  to  friction  can  be  found. 


By  multiplying  the  percent  of  drag  due  to  friction  and  the  percent  increase  in  the  coefficient  of 
friction  due  to  surface  roughness,  the  percent  increase  in  the  total  drag  coefficient  due  to 
roughness  can  be  found. 
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6.  Discussion  of  Tow  Tank  Results 

The  frictional  drag  coefficient  for  the  seven  surfaces  tested  is  shown  on  Figure  12.  For 

fully  rough  surfaces,  Cf  is  constant  for  the  range  of  Reynolds  numbers  tested.  The  smooth 
surface  shows  a  decrease  in  Cf  with  increasing  Reynolds  number.  This  graph  shows  that  all  of 
the  surfaces  tested  are  fully  rough  except  for  possibly  the  80  grit  sandpaper  which  shows  a  slight 
decrease  in  Cf  with  increasing  Rei.  This  does  not  affect  our  results  in  any  way  since  Granville’s 
similarity  theory  can  be  used  for  surfaces  that  do  not  fall  in  the  fully  rough  regime. 
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Using  equation  5,  the  roughness  funetions  for  all  the  surfaees  tested  were  ealeulated.  As 
expeeted,  the  rougher  a  surface  was,  the  larger  its  roughness  function  was.  This  means  that  the 
surface  roughness  with  the  larger  roughness  function  should  also  have  the  highest  drag.  Table  1 
shows  the  roughness  functions  for  the  present  experimental  matrix  used  in  this  investigation. 


Table  1  -  The  roughness  functions  for  the  tow  tank  tests 


Roughness  Functions  (AU*) 

Speed  (m/s) 

80  Grit 

40  grit 

24  Grit 

Fine  Mesh 

Medium  Mesh 

Course  Mesh 

1.86 

5.86 

8.23 

10.08 

6.47 

7.94 

10.03 

2.13 

6.19 

8.65 

10.49 

6.99 

8.36 

10.48 

2.44 

6.53 

8.94 

10.73 

7.24 

8.63 

10.77 

2.74 

6.86 

9.28 

11.09 

7.58 

8.95 

11.09 

3.01 

7.04 

9.46 

11.31 

7.79 

9.15 

11.30 

3.29 

7.25 

9.72 

11.60 

8.07 

9.39 

11.62 

3.59 

7.43 

9.98 

11.82 

8.26 

9.61 

11.84 

Using  equations  (5)  and  (6),  the  roughness  function  A  and  the  roughness  Reynolds 
number,  k^,  were  found  for  the  mesh  and  sandpaper  surfaces.  This  enables  the  full  scale  Cprough 
to  be  determined  as  shown  in  Figure  11.  For  a  DDG-51  at  15  knots,  65  percent  of  the  total  drag 
force  is  frictional.  Thus,  the  total  increase  in  the  drag  coefficient  due  to  roughness  can  be  found. 
Furthermore,  because  power  is  directly  proportional  to  resistance,  the  percent  increase  in 
effective  horsepower,  EFIP,  is  the  same  as  the  total  increase  in  the  drag  coefficient  due  to 
roughness.  Finally,  it  is  known  that  approximately  2.6  million  dollars  a  year  is  spent  on  fuel  for 
each  DDG-5 1 .  By  multiplying  this  amount  by  the  total  drag  coefficient  due  to  roughness,  the 
amount  of  money  spent  in  drag  due  to  surface  roughness  per  DDG  can  be  calculated.  These 


results  are  summarized  in  Table  2. 
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Table  2  -  Results  for  Full  Scale  DDG-51 


Surface 

Roughness 
Height  (pm) 

%  Increase  in 
Frictional  Drag 

%  Increase  in 
EHP 

Cost  of  fuel  due 
to  roughness  in 
millions  of 
dollars 

New  painted 
Ship  Bottom 
Paint* 

130 

6 

4 

0.1 

Fine  Mesh 

380 

60 

39 

1.0 

80  Grit 
Sandpaper 

690 

51 

33 

0.9 

Medium  Mesh 

1170 

74 

48 

1.3 

Course  Mesh 

1420 

104 

67 

1.8 

40  Grit 
Sandpaper 

1800 

78 

50 

1.3 

24  Grit 
Sandpaper 

3600 

104 

67 

1.8 

Fouled 

Ship  Bottom 
Paint* 

5000 

178 

116 

3.0 

*  -  Schultz  (2004) 


The  roughness  penalty  ranges  from  0.9  million  dollars  for  the  smallest  roughness  tested  to  1.8 
million  dollars  for  the  largest  roughness,  which  is  a  significant  portion  of  the  total  amount  of 
money  spent  for  fuel.  The  ability  to  prediet  and  subsequently  reduee  the  drag  due  to  surfaee 
roughness  would  result  in  signifieant  fuel  cost  savings. 

Even  though  ships  will  never  be  covered  with  a  sandpaper  or  mesh  surfaee,  they  still  have 
hull  roughness.  This  hull  roughness  typically  ranges  from  a  peak  to  trough  roughness  height  of 
150  mierons  for  a  newly  painted  ship  to  1  eentimeter  for  a  fouled  ship  eovered  in  barnaeles 
(Sehultz  2003).  The  sandpaper  surfaees  were  measured  in  previous  experiments  using  a  laser 
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profilometer  (Schultz,  private  communication).  This  allowed  for  the  peak  to  trough  roughness 
height  to  be  found.  The  80  grit,  40  grit,  and  24  grit  sandpaper  surfaces  had  maximum  peak  to 
trough  roughness  heights  of  690  microns,  1800  microns,  and  3600  microns  respectively.  The 
mesh  surfaces  were  measured  using  calipers.  The  wire  diameter  was  measured  along  with  the 
pitch,  which  is  the  distance  of  the  space  between  wires.  The  maximum  roughness  height  was 
taken  to  be  twice  the  diameter  of  the  wire  since  the  wires  overlap  each  other.  The  fine  mesh, 
medium  mesh,  and  coarse  mesh  had  maximum  peak  to  trough  roughness  heights  of  380  microns, 
1170,  and  1420  respectively.  The  coarse,  medium,  and  fine  mesh  pitch  to  diameter  ratios  were 
4.5,  2.7,  and  5.1  respectively.  Although  these  surfaces  are  rougher  than  a  brand  new  ship  out  of 
dry  dock  would  have,  they  are  still  considerably  smoother  than  a  fouled  ship  and  thus  significant. 
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7.  Experimental  Procedures  -  Water  tunnel  Tests 

All  direct  measurements  of  the  mean  velocity  profile  and  turbulence  were  done  using  the 
re-circulating  water  tunnel  in  the  Hydromechanics  Laboratory  at  the  United  States  Naval 
Academy.  The  test  section  of  the  water  tunnel  is  0.4  m  x  0.4  m  in  cross  section  by  1 .54  m  long. 
The  flow  is  produced  by  a  four-b laded  axial  impeller  which  is  powered  by  a  55  kW  motor.  This 
motor  can  vary  the  flow  speed  between  0  -  6.0  m/s.  In  this  research  project,  the  speed  of  the 
flow  ranged  froml.5  m/s  to  3.75  m/s.  The  test  plates  (Figure  13)  were  350  mm  wide,  1.54  m 
long,  and  12  mm  thick.  They  were  covered  on  one  side  with  the  roughness  of  interest  which  was 
held  on  using  epoxy.  The  plates  were  then  secured  to  a  fiat  plate  test  fixture  mounted 
horizontally  in  the  bottom  of  the  tunnel.  A  strip  of  36-grit  sandpaper  was  mounted  on  the 
forward  200  mm  of  the  plate  in  order  to  trip  the  developing  boundary  layer  to  turbulent  flow. 
Turbulence  at  the  onset  produces  thicker  boundary  layers  which  allowed  for  easier  measurements 
using  the  LDV  downstream  and  fixes  the  location  of  transition  giving  more  repeatable  results. 
Measurements  were  made  1.35  m  downstream  of  the  leading  edge  of  the  plate.  Previous  studies 
in  the  same  facility  have  shown  that  the  flow  adjusts  back  to  the  surface  condition  of  the  test 
plate  well  upstream  of  the  measurement  location.  Because  the  system  was  closed,  the  growth  of 
the  boundary  layer  created  a  slight  pressure  favorable  gradient  on  the  test  plate. 
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t  location 
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Velocity  and  turbulence  measurements  were  made  using  a  TSI IFA550  (TSI  Inc.,  St. 
Paul,  Minn.,  USA)  two-component,  fiber-optic  laser  Doppler  velocimeter  system.  The  LDV 
used  a  four-beam  arrangement  and  was  operated  in  backscatter  mode.  It  was  mounted  on  a 
three-axis  traverse  unit  (Velmex  Inc.,  Bloomfield,  N.Y.,  USA).  By  crossing  a  laser  light  from 
two  beams  onto  the  moving  fluid,  the  velocity  of  the  fluid  in  one  direction  can  be  measured.  At 
the  crossing  point  the  beams  produce  a  fringe  pattern  due  to  constructive  and  destructive 
interference.  As  the  particle  passes  through  the  fringes,  light  is  scattered.  The  frequency  of  the 
scattered  light  (Doppler  frequency)  is  a  function  of  the  particle  velocity  and  the  fringe  spacing. 
An  example  of  this  can  be  seen  in  Figure  14  (Dantec  Dynamics  2004).  The  scattered  light  is 
measured  by  the  photo  detector,  which  converts  the  light  signal  to  a  voltage.  The  voltage  is  then 
sent  to  a  signal  processor  and  on  to  a  computer  for  high  speed  data  acquisition  and  subsequent 


data  reduction.  Measured  data  is  recorded  at  a  rate  of  100-500  Hz. 
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Figure  14  -  Schematic  of  how  a  laser  measures  the  velocity  of  a  particle 


For  the  current  experiment,  two  components  of  velocity  were  recorded  simultaneously. 
This  requires  two  beam  pairs  that  cross  at  the  same  location.  Each  beam  pair  uses  a  different 
wavelength  of  light.  The  signal  from  each  pair  is  distinguished  using  filters.  The  flow  was 
seeded  with  small  glass  spheres  (10  pm  diameter)  which  are  neutrally  buoyant.  Velocities  are 
obtained  from  the  test  plate  surface  to  the  edge  of  the  boundary  layer  in  the  free  stream  using  40 
logarithmically  spaced  points.  This  enables  multiple  data  points  to  be  taken  close  to  the  wall 
where  the  velocity  gradient  is  large.  A  total  of  20,000  random  samples  per  location  were 
obtained  in  order  to  ensure  converged  turbulence  statistics.  Since  two  components  of  velocity  (4 
beams)  were  obtained  simultaneously,  if  Doppler  bursts  for  the  velocity  measurements  from  each 
pair  of  beams  did  not  fall  within  a  50  pm  coincidence  window,  the  sample  was  rejected. 
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8.  Experimental  Analysis  of  the  Water  Tunnel  Results 

Experiments  were  eondueted  on  three  test  surfaees;  the  smooth  plate,  the  80  grit 

sandpaper,  and  the  fine  mesh  surfaces.  Due  to  the  loss  of  equipment  caused  by  Hurricane  Isabel, 
experiments  were  not  conducted  on  the  other  four  surfaces.  Both  the  smooth  and  sandpaper 
surfaces  were  tested  using  four  different  speeds;  while,  the  fine  mesh  surface  was  only  tested  at 
one  speed.  This  is  shown  in  Table  3. 


Table  3  -  LDV  Test  Matrix 


Reynolds  number  based  on  downstream  distance 
and  freestream  velocity 

Surface 

1.4x10® 

3.4x10® 

4.1x10® 

4.7x10® 

Smooth 

X 

X 

X 

X 

Mesh  -  fine 

X 

X 

X 

X 

Mesh  -  medium 

X 

X 

X 

X 

Mesh  -  coarse 

X 

X 

X 

X 

Sandpaper  -  fine 

X 

X 

X 

X 

Sandpaper  -  medium 

X 

X 

X 

X 

Sandpaper  -  coarse 

X 

X 

X 

X 

X  =  boundary  layer  profile 

Uncertainty  estimates  based  on  multiple  profiles  on  one  surface 


The  turbulent  boundary  layer  velocity  profile  for  the  highest  Reynolds  number  can  be 
seen  in  Figure  15.  This  profile  ofy  vs.  U contains  the  smooth  surface,  the  80  grit  sandpaper 
surface,  and  the  fine  mesh  surface.  As  expected,  the  smooth  surface  profile  is  fuller  than  the 
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rough  surfaces  since  the  smooth  profile  reaches  the  free  stream  velocity  faster.  This  means  that 
that  there  is  a  greater  loss  of  momentum  for  the  rough  profiles,  indicative  of  increased  drag. 


U  (m/s) 

Figure  15  -  Turbulent  boundary  layer  profiles  of  all  three  surfaces  tested 


Measurements  of  boundary  layer  growth  are  the  boundary  layer  thickness,  S, 
displacement  thickness,  S*,  and  the  momentum  thickness,  0.  The  boundary  layer  thickness,  S,  is 
defined  as  the  distance  from  the  surface  where  the  fluid  velocity  is  99.5%  of  the  freestream 
velocity  Ue.  The  displacement  thickness  is  defined  as  the  displacement  of  freestream  fluid  due  to 
boundary  layer  growth  on  the  surface.  The  momentum  thickness  is  a  measurement  of  the  lost 
momentum  due  to  the  surface  drag.  These  parameters  serve  as  measures  of  the  effect  of  surface 
conditions.  The  displacement  thickness,  S*,  is  defined  by  equation  1 1 
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while  momentum  thiekness,  6,  is  defined  by  equation  12. 


W 


u„ 


(12) 


Turbulent  boundary  layer  results  are  often  non-dimensionalized  using  the  frietion 
veloeity,  Ur.  This  is  eonsidered  an  inner  layer  sealing  as  eompared  to  the  free  stream  veloeity 
whieh  is  an  outer  layer  sealing.  In  order  to  find  the  frietion  veloeity,  U/Ug  is  plotted  versus 
\n{yUe  /v)  where  y=yraw+s  and  v  is  the  kinematic  viscosity,  is  the  initial  visual  position  of 
the  laser  at  the  wall,  and  ^  is  a  correction  that  is  added  to  yraw  in  order  to  compensate  for  the  fact 
that  the  laser  was  not  initially  positioned  exactly  at  the  wall  at  the  start  since  the  velocity  at  the 
wall  is  zero.  Furthermore,  with  rough  surfaces,  the  origin  is  not  known  because  it  could  be  at  the 
top  of  the  roughness  peaks,  the  bottom  of  the  roughness  troughs,  or  somewhere  in  between 
(Perry  and  Li  1990).  In  order  to  induce  the  log  law  portion  of  the  inner  region  to  be  linear  in  log 
normal  coordinates,  slight  adjustments  are  made  to  s  and  d (Lewthwaite  et  al.  1985).  This  is 
checked  by  maximizing  R  ,  the  coefficient  of  determination.  An  R  value  can  be  calculated  from 
this  linear  regression  line  which  represents  a  portion  of  the  total  variance  of  the  data  points 
explained  by  the  linear  regression  line.  If  equaled  1,  all  of  the  data  points  would  lie  on  the 
linear  regression  line.  However,  since  s  at  the  moment  is  still  zero,  the  R  value  will  not  be  as 
high  as  it  will  once  a  is  correctly  adjusted.  Once  the  slope  of  the  log  law  region  is  known,  the 
skin  friction  coefficient,  Q,  can  be  found,  as  shown  in  equation  13  (Lewthwaite  et  al.  1985). 
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Cf  =  2k^ 


K^eJ 


& 

V  J 


(13) 


Figure  16  shows  the  skin  frietion  eoeffieient  versus  momentum  thiekness  Reynolds 
number  for  all  surfaces  tested.  The  smooth  wall  skin  friction  coefficient  decreases  with 
increasing  momentum  thickness  Reynolds  number,  whereas  Q remains  nearly  constant  for  the  80 
grit  sandpaper.  As  expected,  both  rough  surfaces  have  higher  C/ values  than  the  smooth  wall.  It 
is  interesting  to  note  that  while  the  mesh  roughness  height  was  smaller  than  the  80  grit 
sandpaper,  it  produces  higher  skin  friction.  Therefore,  skin  friction  is  not  solely  dependent  on 
roughness  height. 


Re^ 


Figure  16  -  Cf  vs,  Ree  (overall  uncertainty  at  95%  confidence  is  4%  for  smooth  surface 

and  7%  for  rough  surfaces) 
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Now  that  Cf,  is  known,  t7ris  calculated  using  equation  14 


U 


T 


where  Ue  is  the  free  stream  veloeity. 


(14) 


The  roughness  funetion,  AU^,  is  found  direetly  from  the  measured  veloeity  profile. 

Sinee  the  y-intereept  for  a  smooth  wall  is  5.0,  the  roughness  funetion  of  a  rough  wall  ean  be 
determined  by  subtraeting  the  known  smooth  wall  intereept  from  the  rough  wall  intereept 
(Lewthwaite  et  al.  1985).  This  ean  be  seen  in  equation  15 

=-lny^  +5.0-AU^  (15) 

K 

whieh  shows  the  mean  veloeity  profile  equation  for  a  rough  wall. 


An  example  of  the  output  of  the  spreadsheet  used  to  perform  all  the  steps  neeessary  to  find  the 
roughness  funetion  ean  be  seen  in  Table  4. 
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Table  4  -  Example  spreadsheet  used  to  find  the  roughness  function 


downstream  position  (mm) 

initial  e  (mm) 

step  size  (mm) 

1350.00 

0.23 

0.01 

lower  limit  constant 

V  (m^/s) 

upper  limit  constant 

1.50 

9.35E-07 

0.15 

CALCULATED  B.L.  PARAMETERS 

Ue  (m/s) 

8  (mm) 

Log  Law  Slope 

3.77 

26.40 

2.44 

Cf 

8*  (mm) 

Log  Law  Intercept 

0.0035 

3.20 

5.00 

u*  (m/s) 

0(mm) 

Log  Law  R 

0.14 

2.53 

0.99995 

roughness  function,  Au+ 

y  lower  limit  (mm) 

y  upper  limit  (mm) 

0.00 

0.37 

3.96 

Table  5  summarizes  all  the  results  of  the  boundary  layer  thiekness,  displaeement 
thickness,  and  momentum  thickness.  As  expected,  there  is  a  significant  increase  in  ^  and  6*  for 


the  rough  surfaces. 
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Table  5  -  Results  from  water  tunnel  tests  (overall  uncertainty  of  roughness  functions  at 

95%  confidence  is  9%) 


3.75  m/s 

8 

(mm) 

8* 

(mm) 

0(mm) 

roughness  function,  AU* 

Smooth  Surface 

26.4 

3.2 

2.53 

0.00 

80  Grit  Sandpaper 

31.9 

5.05 

3.64 

7.72 

Fine  Mesh 

31.1 

5.04 

3.66 

9.00 

3.00  m/s 

8 

(mm) 

8* 

(mm) 

e(mm) 

roughness  function,  AU* 

Smooth  Surface 

26.40 

3.31 

2.60 

0.00 

80  Grit  Sandpaper 

31.40 

5.08 

3.66 

7.22 

Fine  Mesh 

N/A 

N/A 

N/A 

N/A 

2.25  m/s 

8 

(mm) 

8* 

(mm) 

e(mm) 

roughness  function,  AU* 

Smooth  Surface 

27.50 

3.42 

2.64 

0.00 

80  Grit  Sandpaper 

30.80 

4.82 

3.51 

6.43 

Fine  Mesh 

N/A 

N/A 

N/A 

N/A 

1 .5  m/s 

8 

(mm) 

8* 

(mm) 

0(mm) 

roughness  function,  AU* 

Smooth  Surface 

26.00 

3.42 

2.64 

0.00 

80  Grit  Sandpaper 

30.20 

4.98 

3.62 

5.76 

Fine  Mesh 

N/A 

N/A 

N/A 

N/A 
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9.  Discussion  of  Water  Tunnel  Results 

The  mean  veloeity  profiles  at  the  highest  Reynolds  number  can  be  seen  in  Figure  17. 

This  shows  how  the  rough  surface  profiles  of  80  grit  sandpaper  and  fine  mesh  are  shifted 
downward  from  the  smooth  profile;  yet,  they  still  maintain  the  same  shape  as  the  smooth  profile. 
By  using  equation  16,  the  roughness  function  for  the  80  grit  sandpaper  surface  and  the  fine  mesh 
surface  are  determined,  as  listed  in  Table  5.  This  figure  indicates  that  the  mesh  and  sandpaper 
surfaces  tested  display  similar  shifts  from  the  log  law  even  though  the  roughness  type  and  size 
are  different. 


Figure  17  -  Mean  velocity  profiles  for  all  three  surfaces 
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U,-U 

— ^ —  y 

The  velocity  profile,  versus  — ,  is  shown  in  Figure  18.  This  is  called  the 

5 

velocity  defect.  Granville  developed  his  similarity  scaling  on  the  assumption  that  the  velocity 
defect  profile  would  collapse  both  rough  and  smooth  surfaces  alike  (Granville  1987). 
Furthermore,  if  the  roughness  functions  of  the  rough  surfaces  measured  directly  do  not  equal 
those  calculated  from  the  tow  tank  tests,  Granville’s  Similarity  Theory  does  not  apply.  However, 
all  three  surfaces  collapse  on  a  single  curve,  indicating  boundary  layer  similarity. 


0.0  0.2  0.4  0.6  0.8  1.0  1.2  1.4 


(y+£)/S 

Figure  18  -  Velocity  defect  profile 


While  numerous  researchers  have  also  observed  this  collapse  including  Clauser  1954, 
Hama  1954,  Acharya  et  al.  1986,  Granville  1987,  Schultz  1998,  Schultz  &  Myers  2003,  a  few 
studies,  Krogstad  et  al.  1992  and  Krogstad  &Antonia  2001,  noted  differences  in  the  defect 
profiles  for  the  mesh  type  roughness. 
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Finally  a  graph  comparing  the  roughness  type  to  the  roughness  funetion  ean  be  made. 
This  graph  of  Alf  vs  ean  be  seen  in  Figure  19  for  all  of  surfaees  tested  in  both  the  Tow  Tank 
tests  and  the  LDV  experiments,  is  ealled  the  roughness  Reynolds  number  and  represents  the 
non-dimensional  roughness  height  of  a  surfaee.  By  dividing  k,  the  roughness  length  seale  by  the 
viseous  length  seale,  the  roughness  Reynolds  number  is  determined.  The  main  diffieulty  is 
determining  the  proper  roughness  length  seale  to  eollapse  results  from  disparate  surfaees. 


Figure  19  -  Alf  vs  (overall  uncertainty  at  95%  confidence  is  9%) 


The  sandpaper  surfaees  all  lie  on  the  same  line  right  below  the  uniform  sand  eurve.  The 
uniform  sand  eurve  was  first  ereated  by  Nikuradse  who  eondueted  experiments  on  the  inside  of 
eireular  pipes  that  had  tightly  paeked,  uniform  sand  grains  glued  to  the  walls  (Schliehting  1955). 
Uniform  sand  is  defined  as  mono-dispersed  with  uniform  grain  size;  while,  sandpaper  is  poly- 
dispersed  with  a  wider  range  of  grain  sizes.  Sinee  the  sandpaper  surfaees  all  lie  in  a  line  with  the 
same  slope,  a  simple  multiplieation  ean  be  used  to  eollapse  the  sandpaper  surfaees  to  the  uniform 
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sand  curve.  In  order  to  get  the  proper  values  of  k,  the  maximum  peak  to  trough  roughness 
height,  Rt,  was  multiplied  by  a  faetor  of  0.75,  k  =  0.75  Rt  (Sehultz  and  Flack  2003).  Since  Rt  is  a 
measure  of  the  maximum  sandpaper  grain  height,  a  multiplication  factor  must  be  added  to 
aecount  for  the  portion  of  the  surfaee  eovered  with  a  smaller  grain  size.  This  multiplieation 
faetor  cannot  be  determined  analytieally  or  from  a  physieal  measurement  of  the  surfaee,  but  must 
be  determined  experimentally.  Figure  20  shows  how  multiplying  Rt  by  0.75  eollapses  the 
sandpaper  surfaces  onto  the  uniform  sand  curve. 


Figure  20  -  Alt  vs  using  the  sandpaper  scaling  parameter,  A=0,75if, 

Figure  20  also  shows  that  the  tow  tank  tests  and  the  limited  LDV  measurements  eompare 
well.  This  can  clearly  be  seen  with  the  80  grit  sandpaper  surface.  The  LDV  test  results  lie  right 
on  top  of  the  tow  tank  tests  which  is  exactly  what  should  be  happening  if  similarity  laws  hold 
true.  Whereas  a  simple  multiplication  factor  can  be  used  to  eollapse  sandpaper  surfaees;  a 


different  method  must  be  used  with  the  mesh  surfaces. 
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Results  for  the  fine  mesh  and  eourse  mesh  surfaees  do  not  lie  on  the  uniform  sand  eurve, 
however  these  surfaees  follow  a  line  with  a  similar  slope.  For  the  same  value  of  k,  these  surfaees 
have  a  greater  roughness  funetion  than  a  uniform  sand  surfaee.  The  differenees  observed  on  this 
plot  ean  be  explained  by  the  piteh  to  diameter  ratios  and  the  peak  to  trough  roughness  height  of 
the  mesh  surfaees.  The  eoarse  and  fine  mesh  surfaees  have  similar  piteh  to  diameter  ratios, 
whereas  the  medium  mesh  surfaee  has  a  signifieantly  smaller  ratio.  It  seems  that  the  medium 
mesh  surfaee  affeets  the  frietional  drag  in  a  similar  manner  as  uniform  sand  surfaee,  whereas  the 
surfaees  with  larger  piteh  to  diameter  ratios  do  not.  By  performing  a  regression  analysis  of  the 
data,  an  equation  resulted  by  relating  the  roughness  length  seale,  k,  to  the  max  peak  to  trough 
height,  Rt,  and  the  piteh  to  diameter  ratio,  p/d,  as  shown  in  equation  16. 

^-0.2ol  (16) 

d  ) 

Beeause  of  the  limited  range  of  piteh  to  diameter  ratios,  the  funetional  relationship  should  not  be 
expeeted  to  deseribe  all  mesh  surfaees.  Furthermore,  the  graph  whieh  used  this  equation  ean  be 
seen  in  Figure  21 .  Higher  order  polynomial  fits  did  not  yield  signifieantly  better  results 
eonsidering  the  limited  number  of  mesh  surfaees  tested  and  the  uneertainty  in  the  measurement 


data. 
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Figure  21  -  Alf  vs  using  equation  16 


Recently  there  has  been  disagreement  as  to  the  effect  of  surface  geometry  on  the  flow  in 
the  outer  portion  of  the  turbulent  boundary  layer  (Krogstad  and  Antonia  1999,  Antonia  and 
Krogstad  2001).  A  closer  look  at  this  project’s  turbulence  results  will  now  be  made.  There  are 
three  Reynolds  stresses  that  need  to  be  examined.  The  Reynolds  stresses  account  for  the 
additional  stress  due  to  the  mixing  caused  by  turbulent  flow.  The  first  one  measures  the 
Reynolds  stress  parallel  to  the  flow  and  is  also  called  the  axial  stress,  ,  the  second  measures 
the  Reynolds  stress  normal  to  the  flow  and  is  also  called  wall-normal  stress  ,  ,  and  the  third 

measures  the  Reynolds  shear  stress,  -  u'  v'^ .  The  turbulence  stresses,  normalized  by  the 
frictional  velocity,  should  all  collapse  for  y  >  5A:  if  the  surface  roughness  does  not  affect  the 


outer  layer  of  the  boundary  layer.  Figures  22,  23,  and  24  show  the  Reynolds  stresses  start  at  a 
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maximum,  and  then  decrease  to  a  value  of  zero  outside  the  boundary  layer.  This  is  due  to  the 
fact  that  the  largest  changes  in  the  turbulence  occur  near  the  wall. 
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Figure  22  -  Axial  Reynolds  stress 
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Within  the  inner  layer  {yld<  0.2),  the  roughness  type  is  influeneing  the  near  surfaee 
turbulenee.  However,  for  yid  >  0.2,  all  figures  show  excellent  collapse  of  the  Reynolds  stresses. 
This  indicates  that  surface  roughness  does  not  affect  the  outer  layer  of  the  turbulent  boundary 
layer  regardless  of  the  geometry  of  the  roughness. 

The  importance  of  this  result  is  that  current  models  assume  that  the  turbulence  in  the  outer 
layer  is  independent  of  the  wall  roughness.  These  models,  first  published  by  Townsend  in  1976 
are  called  the  Reynolds  number  similarity  hypothesis.  For  the  surfaces  tested,  Townsend’s 
model  holds  true.  It  is  interesting  that  the  wall  can  change  so  much,  from  sandpaper  to  wire 
mesh,  and  the  turbulence  in  the  outer  layer  is  not  affected  when  scaled  by  the  frictional  velocity, 
Ur.  In  addition,  these  results  can  help  researchers  who  study  turbulence  in  the  outer  layer 
because  they  can  now  confidently  disregard  surface  roughness  at  high  Reynolds  numbers. 
Finally,  a  simple  two  layer  empirical  model  based  on  smooth  walls  can  now  be  modified  to 
rough  walls  by  including  skin  friction  to  the  wall  layer. 
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10.  Conclusion 

Seven  different  surfaees  (one  smooth,  three  sandpaper,  and  three  mesh)  were  tested  at 
various  speeds  in  the  115  meter  tow  tank  loeated  in  the  Hydromeehanies  Laboratory  at  the 
United  States  Naval  Aeademy,  and  the  drag  of  eaeh  was  measured.  Using  a  laser  Doppler 
veloeimeter  (LDV),  veloeity  profiles  for  three  surfaees  (one  smooth,  one  sandpaper,  and  one 
mesh)  were  also  measured  in  a  reeireulating  water  tunnel.  The  projeet  produeed  the  following 
results: 

•  The  roughness  funetions  for  the  surfaees  tested  in  the  tow  tank  were  ealeulated  and  the 
full  seale  frietional  drag  eoeffieient  was  predieted  for  a  DDG-5 1  eovered  with  a  similar 
surfaee  roughness.  The  results  indieate  a  signifieant  inerease  in  fuel  eonsumption  with 
inereasing  roughness. 

•  The  roughness  funetions  for  the  veloeity  profiles  measured  with  the  LDV  showed 
exeellent  agreement  with  the  results  from  the  tow  tank  test,  indieating  the  validity  of 
boundary  layer  similarity  laws. 

•  Appropriate  sealing  parameters  for  the  rough  surfaees  were  developed.  The  sandpaper 
sealing  parameter  is  direetly  proportional  to  the  peak  to  trough  roughness  height,  Rt, 
extending  the  results  of  previous  researeh  to  a  larger  range  of  roughness  height.  A  new 
sealing  parameter  for  the  mesh  surfaees  was  developed  whieh  relates  the  piteh  to 
diameter  ratio  as  well  as  the  peak  to  trough  roughness  height  to  the  roughness  funetion. 

•  The  mean  veloeity  profiles  for  the  smooth  and  rough  surfaees  eollapsed  in  defeet 
eoordinates  indieating  similarity  in  the  outer  region  of  the  boundary  layer.  The  Reynolds 
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stresses  also  showed  exeellent  agreement  in  the  outer  layer,  giving  support  to 

Townsend’s  boundary  layer  similarity  hypothesis. 

In  the  future,  additional  mesh  surfaces  with  a  wider  variety  of  pitch  to  diameter  ratios 
should  be  tested  to  enable  the  scaling  parameter  to  be  verified.  Furthermore,  rougher  surfaces 
need  to  be  tested  in  order  to  determine  when  roughness  becomes  too  large  and  starts  affecting  the 
outer  part  of  the  turbulent  boundary  layer.  At  this  point,  current  boundary  layer  models  would 
fail  and  the  similarity  hypothesis  would  no  longer  be  valid. 
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